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The increasing precision of cosmological data provides us with an opportunity to test general 
relativity (GR) on the largest accessible scales. Parameterizing modified gravity models facilitates 
the systematic testing of the predictions of GR, and gives a framework for detecting possible de- 
viations from it. Several different parameterizations have already been suggested, some linked to 
classifications of theories, and others more empirically motivated. Here we describe a particular new 
approach which casts modifications to gravity through two free functions of time and scale, which 
are directly linked to the field equations, but also easy to confront with observational data. We 
compare our approach with other existing methods of parameterizing modied gravity, specifically 
the parameterized post-Friedmann approach and the older method using the parameter set {fi, 7}. 

We explain the connection between our parameters and the physics that is most important for 
generating cosmic microwave background anisotropies. Some qualitative features of this new pa- 
rameterization, and therefore modifications to the gravitational equations of motion, are illustrated 
in a toy model, where the two functions are simply assumed to be constant parameters. 



1. INTRODUCTION 

General Relativity (GR) has been confronted with 
many theoretical and experimental tests since its birth in 
1915. From the gravitational lcnsing experiments in 1919 
El up to the extensive studies and tests of the theory in 
the 1960s and 1970s @-@|, the theory has been confirmed 
observationally and theoretically bolstered in many dif- 
ferent respects. Weak field gravity on experimentally ac- 
cessible scales has been so well tested that there are only 
two remaining directions in which we might find modifi- 
cations to GR: strong field gravity which may be probed 
by studying black holes; and gravity at large scales and 
early times, which is the cosmological arena. 

Cosmology has challenged GR with two, yet to be 
fully understood discoveries: dark matter and dark en- 
ergy 0, § . Along with these two phenomena, the lack of 
renormalizability in GR 0] and the apparently exponen- 
tial expansion in the very early Universe E3] are usually 
taken as signs for the incompleteness of the theory at high 
energies. Due to these shortcomings in GR the study of 
modified gravity has become a broad field. Scalar-tensor 
theories [p, El, f(R) theory El, Horava-Lifshitz 
theory [l5|. multidimensional theories of gravity El, E3 , 
and many other suggestions have been made in the hope 
of finding, or at least deriving some hints for, a fully 
consistent theory that can successfully explain the obser- 
vations and satisfy the theoretical expectations. (Ref. @ 
has an extensive review). 

The new data coming from various experiments such 
as the WMAP and Planck satellite measurements of the 
cosmic microwave background (CMB) anisotropies EH, 
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and the WiggleZ [19( or Baryon Oscillation Spectroscopic 
Survey [2(| measurements of the matter power spectrum, 
provide us with an opportunity to test specific modified 
theories of gravity. However, since there are many differ- 
ent modified theories, all with their own sets of parame- 
ters, there has recently been some effort to come up with 
a way to describe generic modified theories using only 
a few parameters, and to try to constrain those param- 
eters with general theoretical arguments and by direct 
comparison with cosmological data. 

Parametrizing modified theories of gravity with a small 
number of parameters has the benefit of tracking the ef- 
fects of modified gravity on a number of different observ- 
ables consistently and systematically, rather than consid- 
ering the consistency of each single observable with GR 
predictions. This can, at least in principle, lead to the 
acceptance, rejection, or constraining the theory space of 
modified gravity models. 

The parametrized post-Friedmann (PPF) approach, as 
described in Ref. [2l| . is an effective way to parameter- 
ize many of the modified theories of gravity. However, it 
is not really feasible to constrain its more or less dozen 
additional parameters even with the power of Markov 
Chain codes such as CosmoMC (22|; there are just too 
many degrees of freedom to provide useful constraints in 
the general case. In this paper we will describe a slightly 
different way to parametrize modified theories of gravity 
in which we try to retain only a small number of param- 
eters which we then constrain using WMAP 9-year [23| 
and SPT data H. 

In the next section we will describe the formulation of 
this new parametrization, and will show its connection 
with PPF and other approaches in Sec. [3] In Sec. 0] 
we will discuss the results of a numerical analysis using 
CAMB [li] and CosmoMC, and we will conclude the 
paper in Sec. [5] with a brief discussion. 
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2. MODIFIED GRAVITY FORMULATION 

There are two common strategies for modifying grav- 
ity. One can start from the point of view of the La- 
grangian or from the equations of motion. The La- 
grangian seems like the more obvious path for writing 
down specific new theories, where one imagines retain- 
ing some desired symmetries while breaking some others. 
However, the equations of motion provide an easier way 
in practice to parametrize a general theory of modified 
gravity, especially in the case of first order perturbations 
in a cosmological context. 

The evolution of the cosmological background has been 
well tested at different redshift slices, specifically at Big 
Bang nucleosynthesis and at recombination through the 
CMB anisotropics. It therefore seems reasonable to as- 
sume that the background evolution is not affected by 
the gravity modification, with the only background level 
effect being a possible explanation for a fluid behaving 
like dark energy. 

The linearized and modified equations of motion for 
gravity can be written in the following form in a covariant 
theory: 



SG^ = 8irG6T^ + 6U^. 



(1) 



Here, 5G^ V is the perturbed Einstein tensor around a 
background metric, 8T^ V is the first order perturbation 
in the energy-momentum tensor and 5U^ V is the modi- 
fication tensor source from any term that is not already 
embedded in GR. 

Since we will be using CAMB and CosmoMC for 
numerical calculations, we will choose the synchronous 
gauge from now on, and focus only on the spin-0 (scalar) 
perturbations. This will make it much more straightfor- 
ward to adapt the relevant perturbed Boltzmann equa- 
tions. The metric in the synchronous gauge is written 
as 

ds 2 = a 2 (r) [~dr 2 + (<% + h^dx 1 dx 3 ] , 

hij = I d 3 ke ii * s {k i k :j h(k,T) + (kikj - ^%)6?7(fc, r)}, 



k — kk , 



(2) 



where k is the wave vector. Putting this metric into Eq. 
Q] results in the following four equations (2(1 ]: 



Here we have used the following definitions: 

ST° = -6p t a 2 8U$ = k 2 A(k, r), 

5T? = (p + p)Vi , a 2 SU^k 2 B(k,T) 1 



5T! = 3SP 



a 2 SUl = k 2 C{k,T), 



D i:j 5T ij = (p + p)S . a 2 D ij 5U 13 = k 2 D(k, t), 

with p and p being the background energy density and 
pressure, respectively, and a dot representing a derivative 
with respect to r. The factors of k are chosen to make 
the modifying parameters, {A, B,C, D}, dimcnsionlcss. 

Dij is defined as kikj — 

In general, Einstein's equations provide six indepen- 
dent equations. For the case of first order perturbations 
in cosmology, two of these six equations are used for the 
two spin-2 (tensor) degrees of freedom, two of the equa- 
tions are used for the spin-1 (vector) variables and only 
two independent equations are left for the spin-0 (scalar) 
degrees of freedom. This means that Eqs. [3] to |6] are not 
independent and one has to impose two consistency rela- 
tions on this set of four equations. These consistency re- 
lations are of course the energy-momentum conservation 
equation, V ^(T^ + U^) = 0. Assuming that energy con- 
servation holds independently for the conventional fluids, 
V^T^ = 0, (see ref [21| for the strengths and weaknesses 
of such an assumption) one then obtains the following 
two consistency equations: 



2A 
6B + kC 



2A 



2kB 
U 



C 



0: 



YlUB-kD = 0. 



(8) 
(9) 



Here we have defined H = a/a and dropped the argu- 
ments of the functions A to D. 

Eqs. E]to|ni together with Eqs. [5]and[5J show that two 
general functions of space and time would be enough to 
parametrize a wide range of modified theories of gravity. 
This method of course does not provide a test for any 
specific modified theory of gravity. However, given the 
current prejudice that GR is the true theory of gravity 
at low energies (eg, see Ref. [27j for a discussion), the 
main question is whether or not cosmological data can 
distinguish between GR and any other generic theory of 
modified gravity. Clearly, if we found evidence for devi- 
ations from GR, then we would have a parametric way 
of constraining the space of allowed models, and hence 
hone in on the correct theory. 



k 2 V- l;-h= -4i:Ga 2 Sp + k 2 A(k,r); (3) 
2 a 

k<j = 4:irGa 2 (p + p)V + k 2 B(k,r); (4) 

h + 2 - h - 2 k 2 r) = -24 n G a 2 {SP) + k 2 C(k, r);(5) 
a 

h + 6r? + 2 - ( h + Qfj) - 2 k 2 r) = -24ttG a 2 (p + p) S 



+k 2 D(k,r) 



(6) 



3. CONNECTION WITH OTHER METHODS 
OF PARAMETRIZATION 

The parameterization described in the previous sec- 
tion has a very close connection in practice with the 
PPF method explained in Ref. [2lj | . The most impor- 
tant differences arc that we have grouped a number of 
different parameters into a single parameter, and have 
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used the synchronous gauge in Eqs. [3]to[6j The parame- 
ters defined in Sec. [5] are related to the PPF parameters 
according to: 



A(k,r) = A $ + F r + a Q x + -jr-X, 



B(k, T ) = s $ + / r + M + y£ 



(10) 

(ii) 



C(k,r) = C $ + ^l + Jot + yf + 7 oX + yi 



72 

A: 2 ' 



(12) 



51 
~k 2 



(13) 



While the authors of Ref. [2l| have insisted on the modi- 
fications being gauge invariant, it is good to keep in mind 
that there is nothing special about the use of gauge in- 
variant parameters, as is shown in Ref. (28|. The im- 
portant issue is to track the degrees of freedom in the 
equations. There are originally four free functions for 
the spin-0 degrees of freedom in the metric, but the gauge 
freedom can be used to set two of them to zero. Using 
only two gauge invariant functions instead of four, means 
that the gauge freedom has been implicitly used some- 
where to omit the redundant variables. 

All 22 of the parameters on the right hand side of Eqs. 
[TO] to [13] are in fact two-dimensional functions of the wave 
number, k, and time. A hat on a function means that 
it is a gauge invariant quantity. \ 1S the gauge invariant 
form of any extra degree of freedom that can appear, for 
example, in a scalar-tensor theory, or in an f(R) theory 
as a result of a number of conformal transformations (see 
section D.2 of Ref. [2l| for further explanation). $ and 
r are related to the synchronous gauge metric perturba- 
tions through: 



* = Jtf ( h + Qf l + U (h + 617)); 



(14) 
(15) 
(16) 



One needs to add more parameters to the right hand side 
of Eqs. QJJ] to [JJJ] if there is more than one extra degree 
of freedom, or if the equations of motion of the theory 
are higher than second order and the theory cannot be 
conformally transformed into a second order theory. The 
reason this many parameters was introduced in Ref. (2l| 
is that there is a direct connection between these parame- 
ters and the Lagrangians of a number of specific theories, 
like the Horava-Lifshitz, scalar-tensor or Einstein Aether 
theories. Therefore, in principle, constraining these pa- 
rameters is equivalent to constraining the theory space 
of those Lagrangians. 

However, there are a number of issues that may encour- 
age one to consider alternatives to the PPF approach for 



parameterizing modifications to gravity. First of all, it is 
practically impossible to to run a Markov chain code for 
22 two-dimensional functions. One can reduce the num- 
ber of free functions to perhaps 15 using Eqs. [5] and GO 
but there is still a lot of freedom in the problem. The sec- 
ond reason is that the whole power of the PPF method 
lies in distinguishing among a number of classically modi- 
fied theories of gravity that are mostly proven to be either 
theoretically inconsistent, like the Horava-Lifshitz theory 
@, or already ruled out observationally, like TeVeS (at 
least for explaining away dark matter) [231. While it is 
certainly important and useful to check the GR predic- 
tions with the new coming data sets, it does not seem 
reasonable at this stage to stick with the motivation of 
any specific theories. For the moment it therefore seems 
prudent to consider an even simpler approach, as we de- 
scribe here. 

There is another popular parametrization in the litera- 
ture described fully in Ref. [301 an d applied to the case of 
the CMB power spectrum in the same article. This sec- 
ond parametrization is best described in the conformal 
Newtonian gauge, via the following metric: 



ds 2 



2 {t)[-(1 + 2 ip)dr 2 + (1 -2(f))5 ij dx i dx j ]. (17) 



The modifying parameters, {p, 7}, are defined as the fol- 
lowing: 



k 2 ip 



-p(k,a)AirGa 2 {pA + 3(p + p)i:}; (18) 
■y(k,a)ip] = p(k,a)12nGa 2 {p + p)J:. (19) 



Here A = Sp + 3^(1 + p/p)V, and all of the matter 
perturbation quantities are in the Newtonian gauge. 

In order to see the connection between this method 
of parametrization and the one described in the previ- 
ous section, through Eqs. [3] to [H one needs to use the 
modified equations of motion in the Newtonian gauge: 

k 2 <p + 3HU + Hip} = -^Ga 2 Sp + fc 2 A N ; (20) 

k 2 U + H^ = 4TrGa 2 (p + p)kV + k 3 B N ; (21) 



ip) = 12TrGa 2 (p + P)E + k 2 Dj, 



(22) 



The parameters {An, -Bn 7 -Dn} are the modifying func- 
tions in the Newtonian gauge. These parameters are re- 
lated to 7 and p via: 



a = 1 — p: , 
= 7-1, 



(23) 
(24) 



4irGa 2 a{pA + 3(p + p)E} = k 2 (A N - 3—B N 

k 

+D N ), (25) 
(3^^12TTGa 2 aY, + k 2 D N , (26) 

where, one can choose between using the functions 
{An, -Bnj Cn, -On}, along with two constraint equations 
similar to the Eqs. [8] and [9] or using the two parameters 
7 and p and trying to remain consistent in the equations 
of motion. 
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It is argued in Rcf. [31( that the {7, p,} parametrization 
is not capable of parameterizing second order theories in 
the case of an unmodified background and no extra fields. 
To show this they use the fact that, in the absence of 
extra fields, all of the greek coefficients in Eqs. [TO] to 
I13[ i.e. {ao, ...,€2}, have to be zero. Furthermore, they 
argue that in the case of second order theories, Fo and Iq 
have to be zero, and therefore the constraint equations [5] 
and [9] show that Jo and Kq are zero as well. After all of 
this, one can see that equation 1221 can be written as the 
following in this special case: 



k 2 ((j) - ip) = 12nGa 2 (p + P)S + k 2 {D 



t*>- (27) 



The authors then show that the absence of a term pro- 
portional to the metric derivative will lead to an incon- 
sistency. This conclusion is valid only if one assumes that 
that j3 in Eq. [26] is a function of background quantities, 
which usually is not the case. Otherwise, one can use Eq. 
[26] to define /3: 



12 nGa 2 a£ + k 2 (D (, 



Hi 

k 







(28) 



leaving no ambiguity or inconsistency. 1 

It is also claimed in Rcf. (32j that the {7,^} 
parametrization becomes ambiguous on large scales, 
while none of these shortcomings apply to the PPF 
method. However, these criticisms do not seem legiti- 
mate, since, as was shown in this section, there is a direct 
connection between {7, and the PPF parameters. 2 
For any given set of functions for the PPF method, one 
can find a corresponding set of functions {7, /i}, using 
Eqs. QJJ to [13] and [26] that will produce the exact same 
result for any observable quantity. One only needs to 
ensure the use of consistent equations while modifying 
gravity through codes such as CAMB. 

Although we believe that there is no ambiguity in 
the {7, /x} parameterization, we also believe that our 
{A, B, C, D} parameterization can be implemented much 
more easily in codes such as CAMB and CosmoMC. Fur- 
thermore, there is a potential problem for the {7, /1} 
parametrization on the small scales that enter the horizon 
during the radiation domination era. The metric pertur- 
bation ip will oscillate around zero a couple of times for 
these scales and that makes the 7 function blind to any 
modification at those instants of time. This also has the 
potential to lead to numerical instabilities. 



1 Note that this might be troublesome if tp goes to zero at some 
moments of time. This happens for the scales that enter the 
horizon during radiation domination. 

2 In particular there is nothing wrong with the {7, /x} parametriza- 
tion on large scales, since ip is certainly always non zero. 



4. NUMERICAL CALCULATION 

In this section we will constrain the parameterization 
described in Sec. [3] using the CMB anisotropy power 
spectra. We will describe the effects of the modifying pa- 
rameters on the CMB and show the results of numerical 
calculations from CAMB and CosmoMC. 



4.1. Effects of A and B on the CMB power spectra 

Before showing numerical results, we first describe 
some of the physical effects of having non- zero values of A 
or B. So far we have not placed any constraints on these 
quantities, which are in general functions of both space 
and time. There are some effects that can be explicitly 
seen from the equations of motion and energy conserva- 
tion. For example, a positive A enhances the pressure 
perturbation and anisotropic stress, while reducing the 
density perturbation. On the other hand, a positive B 
will enhance the momentum perturbation and reduce the 
pressure perturbation and anisotropic stress. 

There are also some other effects that need a little more 
algebra to see, and we now discuss three such effects: 

Neutrino moments: 

The neutrinos' zeroth and second moments, 
{F„o, F„2}, are coupled to the modifying gravity 
terms according to the Boltzmann equations |26j and 
Eq. |3] to El 



2 • 

F vQ = -kF vl - -h; 



(29) 



F v2 = ^kF vl -^kF v3 + —(h + 6fi). (30) 

Here h is modified according to Eq. [31 and the term 
+ 677 is coupled to A and B through equations [3] and [JJ 



h + 67) 



2k 2 r] + 8irGa 2 5 P 



H 



2A-KGa 2 (p + p) 



2 k 2 A 
H 



6kB. 



V 



(31) 



Therefore, modified gravity can enhance the neutrino sec- 
ond moment. 

Photon moments: 

While the same thing is valid for the photons' second 
moment after decoupling, the situation is different during 
the tight coupling regime. The Thomson scattering rate 
is so high in the tight coupling era that this makes the 
second moment insensitive to gravity. In other words, 
the electromagnetic force is so strong that it does not let 
the photons feel gravity. 

ISW effect: 

The integrated Sachs- Wolfe (ISW) [H effect is pro- 
portional to <j> + ip in the Newtonian gauge. In the syn- 
chronous gauge this is equal to 



ISW <x<p + -ip 



h + 6 fj 
2 k 2 



(32) 
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Here 77 is modified according to Eq. [4] and, therefore, a 
subtle change in the function B(k,r), can have a con- 
siderable influence on the ISW effect. Fig. Q] shows the 
effects of a constant non-zero A and B on the ISW effect. 
For the case of a constant non-zero B, the ISW effect is 
always present, since the time derivative of the potential 
is constantly sourced by this function. This will result in 
more power on all scales, including the tail of the CMB 
power spectrum (see Fig. [2]). Fig. [2] shows that if a non- 
zero B is favoured by CMB data, it is mostly due to the 
large i's, (I > 1500), and it is due to its anti-damping 
behaviour. 
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FIG. 1: The effects of a constant, non-zero A and B on the 
ISW effect. This plot shows the ISW effect for a scale of 
k = 0.21Mpc -1 . 




FIG. 2: Effects of a constant B on the CMB power spectra. 
The plot shows the difference in power for the case of zero 
B minus the best fit non-zero B, using WMAP9 and SPT12 
data, while keeping all the rest of the parameters the same. 
The error band plotted is based on the reported error on the 
binned CMB power spectra from the WMAP9 [H|] and SPT12 
[H] groups. 



Fig. [3] shows the CMB power spectra for the case of 
a constant but non-zero A and B. Note how a constant 
non-zero B brings the tail of the spectrum up. One might 
also point to the degeneracy of A and the initial ampli- 
tude mostly by looking at the height of the peaks. 
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FIG. 3: The effects of a constant, non-zero A and B on the 
CMB power spectra. One can see that the two parameters 
have quite different effects. 

Matter overdensity: 

The Boltzmann equation for cold dark matter overden- 
sity in the synchronous gauge reads [26j |: 



^CDM 



Sp 



CDM 



PCDM 



(33) 



Using Eqs. [3j Q] and the Friedmann equation, V. 2 = 

SnGa 2 _ . , . 

p, and assuming a matter dominated Universe 

.3 

with no baryons, one gets the following equation for the 
cold dark matter overdensity: 




k 2 r] 



3HS. 



CDM = 



mathcalH 



(34) 



The above equation clearly shows the role of A and B as 
driving forces for matter overdensity. The fc 3 prefactor 
makes the first term on the right hand side dominant on 
the small scales and will therefore have a significant effect 
on the matter fluctuation amplitude at late times. Mat- 
ter power spectrum will therefore put strong constraints 
on modified gravity models. 



4.2. Markov Chain constraints on A and B 

Since A(k, r) and B(k, r) are free functions, we need to 
choose some simple cases to investigate. We choose here 
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to focus on the simple cases of A and B being separate 
constants (i.e. independent of both scale and time). We 
do not claim that this is in any sense a preferred choice 
- we simply have to pick something tractable. With 
better data one can imagine constraining a larger set of 
parameters, for example describing A and B as piecewise 
constants. 

We used CosmoMC to constrain constant A and B, 
together with the WMAP-9 H and SPT12 [U CMB 
data. The amplitudes of the foregrounds were added as 
additional parameters and were marginalized over for the 
case of SPT12. The resulting constraints and distribu- 
tions are shown in Fig. 2J 

WMAP9 WMAP9+SPT12 



~ 0.0005 I 





constant A 



-0.02 -0.015-0.01 -0.005 0.005 0.01 0.015 0.02 

Constant A 



-0.02-0.015 -0.01 -0.005 0.005 0.01 0.015 1 

Constant A 



FIG. 4: The 68 and 95 percent contours of the constants A 
and B using WMAP9-year data alone (left) and WAMP9 + 
SPT12 (right). 

There is a strong (anti-)correlation between A and 
the initial amplitude of the scalar perturbations, and 
therefore the constraints on A arc quite broad. Two- 
dimensional contour plots of A versus As, the initial am- 
plitude of scalar curvature perturbations, are shown in 
Fig. El 

The mean of the likelihood and 68% confidence interval 
for the six cosmological parameters together with A and 
B are tabulated in Table U Notice that while WMAP9 
is fairly consistent with the conventional ACDM model, 
there is a 3<r preference for non-zero B when one adds 
the SPT12 data. 



FIG. 5: The strong anti-correlation between parameter A and 
the initial amplitude, A s makes the constraints on any one of 
these two parameters weaker. 



if B was an oscillating function of time, changing sign 
from time to time, its total effect on the CMB power 
spectra would become weaker and the constraints would 
be broader. According to Eq. [4] a constant B will change 
r] monotonically, while the effect of an oscillating B will 
partially cancel some of the time. 

By comparing the results of the MCMC analysis for 
WMAP and SPT, it seems like a non-zero B is somewhat 
favoured by data. Essentially what is happening here is 
that a mild tension between WMAP and SPT is resolved 
by adding an appropriate degree of freedom. One can 
see from Fig. [3] that crudely, a constant B changes the 
normalization of the large-scale vs. small-scale power 
spectrum. However, it will take significantly better data 
to determine whether the improved fit genuinely looks 
like a modification of gravity, or has a more mundane 
explanation. In particular, as was mentioned before, the 
amplitudes of the foreground power spectra were consid- 
ered unknown parameters and were fitted and marginal- 
ized over to obtain the distributions for A and B here. 
Since B is slightly correlated with all of these amplitudes 
(and mostly with the point source shot- noise amplitude), 
then upcoming data from the Planck satellite certainly 
will be helpful. 
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5. DISCUSSION 

Since a constant A is essentially degenerate with the 
initial amplitude of the primordial fluctuations, the CMB 
cannot constrain this parameter. On the other hand, B 
seems to be fairly well constrained by the data. However, 
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Parameter WMAP9 WMAP9+SPT12 



100ST2 b /i 2 2.22 ±0.05 2.10 ±0.03 

tt DM h 2 0.118 ±0.005 0.122 ±0.005 

100 1.038 ±0.002 1.040 ±0.001 

t 0.086 ±0.014 0.076 ±0.012 

log{10 w A B ) 3.1 ±0.2 3.1 ±0.2 

n B 0.96 ±0.01 0.93 ±0.01 

100A 1± 10 0.02 ±10 

1000B 0.44 ±0.41 0.74 ±0.24 



TABLE I: The mean likelihood values together with the 68% confidence interval for the usual six cosmological parameters, 
together with constant A and B. 
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